Electronic structure and rebonding in the onion-like As@Nii 2 @As 2 o cluster 
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We present the first ab initio study of the geometry, electronic structure, charged states, bonding 
and vibrational modes of the recently synthesized fullerene-like As@Nii2@As2o cluster which has 
icosahedral point symmetry [Science, 300, 778 (2003)]. We show that the molecule is vibrationally 
stable and will be electronically most stable in its -3 oxidation state in the condensed phase and 
in —2 state in the gas phase. We examine the bonding in this unusually structured molecule from 
charge transfer between atoms, infrared and Raman spectra, and charge density isosurfaces. 
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Since the discovery that carbon can form fullerene-like 
structures of high symmetry 0, Q , there has been con- 
siderable speculation, both experimentally and theoreti- 
cally, as to whether symmetric fullerene-like clusters can 
be created from other elements 0, . Solids comprised 
of such inorganic sub-units can have important applica- 
tions in nanotechnology and materials science. An impor- 
tant step towards this goal was realized with the recent 
synthesis of the [(C4H9)4P +1 ]3[As@Nii2@As2o]~ 3 salt by 
Moses et al. Q where the [As@Nii 2 @As 2 o]" 3 anion (Cf. 
Fig. la) is comprised of an icosahedral (1^) As@Nii2 
cluster with an arsenic atom core (Cf. Fig. lb) encapsu- 
lated by a dodecahedral fullerene-like AS20 cage (Cf. Fig. 
lc) to give an onion-skin-like [As@Nii2@As2o]~ 3 cluster 
with icosahedral point symmetry. As discussed by Moses 
et al. 5], clusters with this structure form the inner 
core of a class of intermetallic compounds which were 
first identified by Bergman et al. @. Such cores might 
actually be responsible for the formation of quasicrys- 
tals. The encapsulating AS20 is shows only structural 
resemblance to a carbon-based fullerene. Electronically, 
the atoms in AS20 are sp 3 -hybridized with three bond- 
ing orbitals along the surface of the dodecahedron and a 
fourth orbital containing a lone pair of electrons pointing 
outward from the cage. In the carbon-based fullcrencs 
the carbon atoms are sp 2 -hybridized forming a strong 
covalently-bonded network of cr-bonds surrounded by a 
sheath of 7r-electrons in bonding orbitals. 

The As@Nii2@As2o cluster is fascinating in that it rep- 
resents another example of what Miiller has called the 
"beauty of symmetry" where the interpenetrating na- 
ture of the AS20 dodecahedron and the As@Nii2 icosa- 
hedron yield a polyhedra of sixty triangular faces 0. 
Further, it demonstrates the first experimental confir- 
mation of a completely non-carbon fullerene-like dodec- 
ahedron, As2o- At this stage, very little is known about 
the molecule and the only theoretical study to date has 
calculated the molecular orbital energies with the semi- 
empirical extended Hiickel method |5|. In this commu- 



nication, we present results of our first-principles density 
functional theory (DFT) studies on the isolated neutral 
and anionic states of the As@Nii2@As2o cluster and its 
sub-units, AS20 and As@Nii2 , to determine their respec- 
tive electronic structure, vibrational stability, and chem- 
ical bonding patterns. In particular, our calculations 
demonstrate that a rehybridization of the the covalent 
As- As bonds present in the isolated AS20 cage takes place 
to form new As-Ni bonds between the inner As@Nii2 
cluster and the outer AS20 cage. 

Our DFT calculations were performed at the all- 
electron level with the generalized gradient approxima- 
tion (GGA) @ to describe the exchange-correlation ef- 
fects. We have performed these calculations with T/j 
symmetry but we find that the optimized geometry has 
lh symmetry which was verified from the positions of the 
atoms. All calculations have been performed using the 
NRLMOL package which employs a Gaussian basis 
set where the exponentials are optimized for each atom. 

The optimized equilibrium geometries of the isolated 
As@Nii2@As2o molecule along with its two sub-units are 
shown in Fig. 1. For the free standing icosahedral 
As@Nii2 , the core As-Ni distance is 2.36 A while the Ni- 
Ni distance is 2.48 A. In the AS20 cage, the As- As distance 
is 2.5 A which is in good agreement with earlier calcu- 
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FIG. 1: Structures of a)As@Nii 2 @As 2 o , b) As@Nii 2 , and 

c) AS20 • 



2 



55 



As(s) p 




1 1 1 1 , i 1 i 1 1 1 h 

As(c) p 

A A 


A 

^-H 1 1 1 1 


As(c) s 


! — 










— 1 — 1 — 1— ^ — 1 — 1 — 1 — 1 

Nid 

i - i a A- W I 


-a 




,/ 


V/ 


1 1 1 1 1 

- , ~v— 


Ni p 

, /-W 






1 1 1 




A Total minority 

'! ajl a/\ i 








\r 





-il -10 -9 -8 -7 -6 -5 -4 -3 -2 



Energy (eV) 

FIG. 2: The total and atom projected partial density of states 
in arbitrary units of As@Nii2@As2o of the minority channel. 
The three upper panel shows the projected density on the 
central (c) and outer shell (s) As atoms. The Ni p DOS is 
magnified for presentation purpose. 

lated results ^(j- Upon formation of the As@Nii2@As2o 
structure, the central As-Ni, Ni-Ni, and the outer cage 
As- As distances increase to 2.61, 2.74, and 2.8 A, respec- 
tively. The Ni to outer layer As distance is 2.43 A. The 
interatomic distances increase by about 10-12 % upon 
formation of the As@Ni 12 @As2o molecule as compared 
to isolated sub-units. These bond lengths are in excel- 
lent agreement with the experimental values (Cf. Table 
HJ. We note that while our calculated bond lengths are 
for the neutral molecule, the experimental values are for 
the molecule in the crystalline environment where it is in 
—3 charged state and the perfect I/, symmetry is weakly 
broken by standard crystal field effects and by the posi- 
tively charged ligands. Our calculations have shown that 
the structural changes induced by the excess charge are 
small - from 0.4 to 1 %. 

The calculated atomization energy of the free-standing 
As 20 unit is 2.79 eV and that of the As@Nii 2 is 3.05 eV. 
The As@Nii2@As2o molecule has atomization energy of 
3.7 eV and is stable by about 26 eV with respect to its 
two sub-units. Its ionization potential is 6.4 eV. 

The calculated total and atom projected density of 
states (DOS) of the isolated neutral molecule are shown 
in Fig. [21 Since the majority states are completely occu- 
pied, we show the DOS for minority states responsible for 
chemical interaction. Both the highest occupied molec- 
ular orbital (HOMO) and the lowest unoccupied molec- 
ular orbital belong to the minority states. The HOMO 
has H„ symmetry, while the LUMO has Ti„ symmetry. 
The eigenvalues of HOMO and LUMO along with the 
Hubbard U values are presented in Table [I] The neutral 
cluster has a spin magnetic moment of 3 fiB- Reduction 



to the —3 charged state quenches the magnetic moment 
and opens up a large HOMO-LUMO gap, as is evident 
from the Fig. [21 This in turn stabilizes the molecule 
by 2.4 eV. The HOMO-LUMO gap in the neutral and 
triply charged state are 0.115 and 1.45 eV, respectively. 
This observation is in accord with the magnetic suscep- 
tibility measurement as well as with the Hiickel model 
prediction |5j. The DOS of the As@Nii2 molecule (not 
shown) shows As s states in the excited levels which be- 
come occupied in the As@Nii2@As2o complex. The iso- 
lated As 2 o is stable with a large HOMO-LUMO gap of 
1.44 eV. In As@Nii 2 @As 2 o , Ni d minority spin states 
are filled and the LUMO have a predominantly p char- 
acter. An analysis of the atom projected DOS and the 
orbital densities shows that both the HOMO and LUMO 
are centered around the Ni atoms. Our calculation on —1 
and —2 charged states show that the isolated molecule is 
most stable in the —2 charged state. Although the elec- 
tronic states are completely filled in the triply charged 
anion, the electrostatic repulsion due to the extra elec- 
trons makes the molecule energetically unfavorable in 
the gas phase. The first and second electron affinities 
of the As@Nii2@As 20 are 3.47 and 0.74 eV respectively. 
The accuracy of the electron affinities is subject to the 
exchange-correlation functional used . The low bind- 
ing energy of the second electron makes it susceptible to 
detachment in the laser desorption experiments. This is 
consistent with the observation that the mass spectrum 
shows peaks for the singly charged anions pj. 

To determine qualitatively the nature of bonding be- 
tween Ni and the As shell, we have integrated the charge 
density inside non-overlapping spheres of arbitrary radii 
around each atom . These values in units of e are given in 
Table [H] and show a very small amount of charge trans- 
fer upon formation of the onion-like molecule especially 
around the outer shell As atoms. The same feature is 
seen in various charged states of the molecule. There 
could be two possible causes for the small charge trans- 
fer between the atoms. First, the interaction between the 
AS20 and As@Nii2 cages may be weak and secondly, the 
environment around the outer As atoms may be similar 
in As 20 and in As@Nii 2 @As 2 o ■ 

The strength of interaction between the two cages 
will also be reflected in the vibrational modes of 
As@Nii2@As2o and its infrared (IR) and Raman activity. 
If the bonding between the cages is weak, then the vibra- 
tional modes of the As@Nii2@As 2 o will be quite similar 
to those of isolated AS20 and As@Nii 2 . 

The vibrational frequencies for As@Nii2@As2o and 
both of its sub-units are real indicating that they are 
vibrationally stable. The vibrational frequencies of this 
massive molecule range between 42 to 305 cm -1 . The 
frequencies of the AS20 and As@Nii2 range from 63 — 248 
cm -1 and from 104 — 305 cm -1 , respectively. The mean 
polarizability of As@Nii2@As 2 o is 25.7 bohr 3 per atom 
which indicates that this molecule may be useful for ap- 
plication in optical devices. From the polarizability to 
molecular volume ratio, we find that the polarizability 
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TABLE I: The experimental and calculated values of bond lengths are presented in A. Note that (c) refers to the central As 

atom. The eigenvalues of the HOMO and the LUMO levels are given in Hartrees. The last column shows the Hubbard U 
parameter in Hartree/electron 2 . 

Method As - As Ni - Ni As - Ni As(c) - Ni Ehomo E LU mo U 
T3xpT 2/75 2^40 2/36 

As 20 2.50 -0.1974 -0.1443 0.1052 

As@Nii2 2.48 2.36 -0.1445 -0.1384 0.1266 

As@Nii2@As 20 2.80 2.74 2.43 2.61 -0.1827 -0.1789 0.1029 
[As@Nii2@As2 r 3 2.80 2.72 2.44 2.59 



of the As@Ni 12 @As 20 is about 10% larger than in C 60 
molecule. The IR absorption (top three panels) and the 
Raman spectra (bottom three panels) of As@Nii2@As2o 
together with its isolated sub-units are presented in Fig. 
El The projection of the spectra of As@Nii2@As2o on 
the As and Ni cages are also shown in the same figure. 
The As@Nii2@As2o has ten optical modes - four IR and 
six Raman active modes. These modes and their corre- 
sponding symmetries are tabulated in Table ITTT1 

The triply degenerate IR active modes belong to the 
Ti„ symmetry. The IR spectrum of As@Nii2@As 2 o is 
dominated by a large peak at 265 cm' 1 with three smaller 
peaks at 84, 163, and 212 cm -1 . The AS20 cage shows 
little IR activity with a single peak at 167 cm -1 which 
disappears upon encapsulation. The low frequency IR 
mode of As@Nii2@As2o at 84 cm' 1 appears to be due 
to the rocking motion of the As@Ni 12 in which the AS20 
cage also gets mildly distorted. This rocking motion is 
also responsible for the peak seen at 143 cm -1 in the IR 
spectra of As@Nii2 . The small peak at 163 cm -1 in the 
IR of As@Nii2@As2o is due to the motion of the central 
As atom. The most dominant peak in the spectrum of 
As@Nii2@As2o at 265 cm' 1 is caused by torsional motion 
mostly of the Ni cage. This peak can not be correlated 
with the peak at 305 cm -1 in the As@Ni 12 spectra. 

The Raman active modes have Ai g and H ff symme- 
try which are one and five-fold degenerate. The lowest 
modes of the AS20 and As@Nii2 couple and form the 
two peaks at 67 cm -1 and 102 cm' 1 upon encapsula- 
tion. The increased intensity of the lowest Raman peak 
is a result of increase in polarizability derivative upon en- 
capsulation. In fact, we find that the ratio of the Raman 
scattering intensities of As@Nii 2 @As2o and AS20 scales 
exactly as the square of the ratio of their radii. Similar 
agreement between the intensities and the ratios of the 
radii is also observed for As@Nii2@As2o and As@Nii2 . 



These modes are 5-fold symmetric involving stretching 
of the molecule diameter. The mode at 142 cm' 1 is a 
breathing mode and correlates to the breathing mode of 
AS20 at 152 cm -1 . The As@Nii2 shows a Raman peak 
due to breathing motion at 263 cm' 1 . This mode and 
the breathing mode of the AS20 cage at 152 cm -1 couple 
together to form a peak at 227 cm' 1 which is basically 
an asynchronous breathing motion with motions of the 
Ni and As cage being in opposite phase. 

The absence of a one-to-one mapping between the vi- 
brational modes of As@Nii2@As2o and its sub-units in 
the IR and Raman spectra indicates that the interaction 
between these sub-units is non-negligible. We do not see 
any pure "rattling" motion of As@Nii2 inside the AS20 
cage. On the contrary, there seems to be a strong cou- 
pling of the vibrational modes between the As@Nii2 and 
the AS20 cages. 

The same feature of strong interaction is also seen 
in the isodensity surfaces of the As@Nii 2 @As2o clus- 
ter. In Fig. 01 we present the isovalued surfaces of 
charge density in the isolated AS20 and in the encapsu- 
lated As@Nii2@As2o cluster. The density surface in AS20 
shows covalent bonding between the As atoms. The As 
atoms in AS20 are sp 3 hybridized and form covalent a- 
bonds with its three nearest As neighbours. The fourth 
orbital contains a lone pair of electrons pointing radially 
outward. However, in the encapsulated form, each of the 
outer As atoms have three nearest neighbor Ni atoms and 
form strong covalent bonds with the Ni d orbitals rather 
than with other As atoms. We refer to this phenomenon 
as rebonding of the As atoms. Thus each As atom has 8 
surrounding valence electrons and each of the Ni atoms 
form such covalent bonds with its five nearest As atoms. 
The Ni atoms occupy a position near the center of the 
pentagonal faces but slightly out of plane and since the 
As-Ni bonds are stronger than Ni-Ni bonds, the Ni cage 
expands. The expansion of the As cage occurs possi- 



TABLE II: The charge within a sphere of radius 2.28 a.u. 

and 2.17 a.u. around the As and Ni atoms, respectively. The TABLE III: The IR and Raman active frequencies of 
symbols (c) and (s) refer to the central and outer shell. As@Nii 2 @As 2 o • The values are in cm" 1 . 
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FIG. 3: The IR and Raman spectra for the As@Nii2@As2o 
molecule and its sub-units (see text for details). 

bly due to the kinetic energy repulsion. It may be men- 
tioned here that the most favorable dissociation channel 
for isolated AS20 is AS4 T(| whereas the laser desorption 
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mass spectroscopy shows peaks for all As@Nil2@As x for 
x=l,20 Since in the encapsulated structure, each As 
is now bonded strongly with three nearest Ni atoms and 
the As-As bonds are much weaker than in AS20 , the 
laser desorption experiment breaks away each As atom 
individually rather than in AS4. 

In conclusion, we find that the isolated As@Nii2@As2o 
cluster has 1^ symmetry and is vibrationally stable. We 
find that the molecule will be electronically stable in the 
—3 charged state in the condensed phase due to the open- 
ing of a large HOMO-LUMO gap. The electrostatic re- 
pulsion makes it more stable in the —2 charged state 
when in the gas phase. Our analysis of the IR and Ra- 
man spectra of the molecule vis-d-vis its sub-units and 

the electronic structure shows that the outer cage As- 
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FIG. 4: The iso-surface of the total charge density in AS20 
and in As@Nii 2 @As 2 o cluster. 

As bonding changes dramatically upon encapsulation of 
As@Nii2 . There exists strong covalent bonding between 
the outer As atoms and inner Ni atoms, in contrast to the 
Huckel-based prediction of Ref. 5. Our study provides 
an understanding of the bonding and will be useful for 
further experimental characterization. 
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